Quantum quench in ID: Coherent inhomogeneity amphfication and 'supersohtons' 
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We study a quantum quench in a ID system possessing Luttinger liquid (LL) and Mott insulating 
ground states before and after the quench, respectively We show that the quench induces power law 
amplification in time of any particle density inhomogeneity in the initial LL ground state. The scaling 
exponent is set by the fractionalization of the LL quasiparticle number relative to the insulator. As 
an illustration, we consider the traveling density waves launched from an initial localized density 
bump. While these waves exhibit a particular rigid shape, their amplitudes grow without bound. 

PACS numbers: 



The shattering of cold glass in hot water is but one of 
many spectacular effects that can be induced by a rapid 
thermal quench in classical media. What happens when 
an isolated quantum phase of matter is subject to a sud- 
den, violent deformation of its system Hamiltonian (a 
'quantum quench')? This question is now under vigor- 
ous investigation in cold atomic gases [H-Q- Long-time, 
out-of-equilibrium physics already observed in gases con- 
fined to one [H, two [3|, and three Q spatial dimensions 
includes oscillatory collapse and revival phenomena [H, Q 
and topological defect formation [sl, [H . 

In this Letter, we study interaction quenches in one- 
dimensional (ID) quantum many body systems. Prior 
theory assuming spatially uniform dynamics has consid- 
ered the post-quench distribution of quasiparticles [6|, 
correlation functions 0, [H , thermahzation [gI, [qI , quan- 
tum critical scaling jlO|, etc. On the other hand, the 
stability of homogeneous solutions with respect to the 
spontaneous eruption of spatial non-uniformity is by no 
means guaranteed, due to the coupling between modes 
with different momenta and the extensive quantity of en- 
ergy injected into the system by the quench. Indeed, 
homogeneous external perturbations are known to gen- 
erate large spatial modulations in a variety of physical 
contexts [5|, [uj . We show here that quantum quenches 
can produce strongly inhomogeneous states via a mech- 
anism that is ubiquitous in ID. 

We consider quenches across a quantum critical point, 
with initial (pre-) and final (post-quench) Hamiltoni- 
ans possessing Luttinger liquid (LL) and Mott insulator 
ground states, respectively. Specifically, we quench into 
the insulating phase of the quantum sine Gordon model 
at the "Luther-Emery" (LE) point [1, [IqI, [iMi , where 
we are able to determine the dynamics analytically. The 
pre-quench ground state has an inhomogeneous density 
profile po{x)^ which acts as a "seed" generating fluctu- 
ations in the space-time dynamics of local observables 
[Til . [l6| . We find that an arbitrarily small deviation of 



po{x) from a constant is dynamically amplified by the 
time evolution, see e.g. Figs. [1] and [2l We argue that 
the mechanism responsible for the amplification is quasi- 
particle fractionalization^ a generic attribute of gapless 



interacting particles in ID [12|, ll7|. We further illustrate 



the amplification effect for a localized (Gaussian) initial 
density "bump." This bump gives rise to a pair of non- 
dispersive, non-interacting density waves that exhibit a 
rigid shape, with amplitudes that grow in time as a power 
law. We have dubbed these traveling density waves 'su- 
persolitons'; an example is depicted in Fig. [2l 

Specifically, for the Fourier transform p(t, k) of the den- 
sity operator expectation value p(t, x), we find the follow- 
ing exact asymptotic result, valid in the long time limit: 

p{t,k) 



Po(^) 



COS(HO - Aa mty/^ COS (\k\t' + ^) , (1) 



where Aa is a non-universal, /c-independent constant and 




FIG. 1: Space-time evolution of the right-moving number 
density pR after Luttinger liquid to Mott insulator quench, 
demonstrating the instability of spatially uniform dynamics; 
fainter (bolder) traces depict earlier (later) times. An in- 
finitesimally small initial density inhomogeneity grows with- 
out bound. The figure is obtained from Eq. ([T]) with a — 0.8, 
Aa = 4.7, and an initial density profile po{x) given by a sum of 
150 cosines with random amplitudes, phases, and wavenum- 
bers. Amplification occurs for any non-zero cr, corresponding 
to a non-zero fractionalization of the initial LL quasiparticles 
with respect to the insulator. 
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t' = t/K^ where K = 1/4 locates the LE point (see 
below); the quench is performed dXt' = 0. The exponent 
cr in Eq. ([1]) is determined by the relative fractionalization 
of the LL quasiparticle number with respect to the Mott 
insulator, 



cr = {K/2K^K/2K) - 1, 



(2) 



where K is the Luttinger parameter characterizing the 
initial Hamiltonian. Eq. ([T]) implies that the density 
splits into non-dispersing left- and right-moving compo- 
nents, p(t, x) = Pr{x — t') pl{x -\- 1'). Interestingly, the 
long time response is linear in pQ and enhanced at shorter 
wavelengths due to the fractional derivative {\k\^^'^) fac- 
tor. For cr > 0, the fluctuations of pr^l are continuously 
amplified by the quench. The effect is demonstrated in 

Fig.m 

In the rest of this Letter, we will explain the setup and 
calculations leading to Eq. ([T]). Before the quench, our 
cold atom system is assumed to reside in the ground state 
|0)^^ of the LL Hamiltonian 



-I 



dx 




Pojx) dO 
dx 



(3) 

where v is the sound velocity, K is the Luttinger pa- 
rameter, and po{x)/q is an external chemical potential, 
with q = K/v7T. The Hamiltonian in Eq. (jSj governs 
the low-energy, long- wavelength physics of many gapless 
ID cold atomic and condensed matter quantum systems 
12, ITsI; in this paper, we have in mind a ID optical 
lattice gas of spin-polarized, neutral Fermi atoms, but 
other interpretations are possible. The short-ranged in- 
teratomic interactions determine v and K; repulsive (at- 
tractive) interactions correspond to < 1 {K > 1), 
while the free Fermi gas has K = 1 and v equal to the 
bare Fermi velocity. The boson fields (j) and encode 
fluctuations of the long wavelength fermion number den- 
sity :p : and current : J : on top of the filled Fermi sea via 
^/7^ :p:= dO/dx and ^/7T :J:= d(j)/dx^ where :. . . : denotes 
normal-ordering with respect to the homogeneous ground 
state |0)^^^Q. These satisfy the commutation relations 

[:p{x) :,: J(x') :] = -{i/7T){d/dx)S{x - x'). The static 
chemical potential in Eq. (|3]) allows us to "write" an ar- 
bitrar 
12 



ary 



density profile into |0) via the axial anomaly 



(0| ■.p{x):\0)^=po{x), 



■J{x):\0)^=0. (4) 



With our system initially prepared in the LL ground 
state |0)^^, we perform the quench at time t = 0. The 
dynamics for t > are generated by the translationally 
invariant, "final state" Hamiltonian Hf, which favors a 
gapped, Mott insulating ground state. Specifically, Hf is 



the Hamiltonian of the quantum sine Gordon model, 

/-\2 /.\2 




cos 



(5) 



In Eq. ([5]) we have expressed Hf in terms of the canon- 
ically rescaled boson variables <l> = ^J^f(t) and 6 = 
0/ ^Kf. The Mott gap-inducing interparticle interac- 
tions set the parameters M and Kf. In the context of a 
Fermi lattice gas at commensurate filling, the "Luttinger 
parameter" Kf characterizes pure forward scattering, 
while M gives the strength of backward scattering Umk- 
lapp interactions; a is a cutoff-dependent length scale. 
The ground state of Hf is gapped for arbitrarily small 
M over the regime < Kf < 1/2, in which the quan- 



tum sine Gordon model is integrable [12|. The solitons 
and antisolitons of the classical sine Gordon equation ap- 
pear as massive Dirac fermions in the quantum version 
[201 . Solitons repel antisolitons for 1/4 < Kf < 1/2 and 
attract them for < Kf < 1/4; in the latter case, ad- 
ditional bosonic bound states {breathers) appear in the 
spectrum. We choose to quench to the boundary between 
these two regimes, where Kf = K = 1/4. At this spe- 
cial "Luther-Emery" point, the interactions between the 
quantum solitons switch off, and Hf can be refermionized 
|14| in terms of a massive non-interacting soliton field ^, 



Hf 



-.1 



d 

dx 



(6) 



In this equation, ^ is a two-component Dirac fermion 
that is related to the boson fields in Eq. ^ via the 
bosonization identity, oc exp[iy^(<l> ± O)]; a^'^ are 

Pauli matrices in the standard basis. The mass gap M 
in Eq. ([6]) is a non-universal, cutoff-dependent quantity. 
It is instructive to rewrite Hi [Eq. (|3])] in terms of ^, 



H, 



J dx ji;^^ 



2K^ 



dx J 2q 
[K'^ -K'^] : :: :| , (7) 



where v = Kv/K. Comparing Eqs. (|6j) and ([7]), we 
see that the quench with = ^ is special. For this 
case only ("non-interacting" quench), the quasiparticles 
of the initial and final Hamiltonians are in one-to-one cor- 
respondence. At any other value of 7^ ^ ("interacting 
quench"), an elementary excitation of the initial state 
carries a fraction of the final state quasiparticle number; 
that is, the "quasiparticle" excitations of the initial LL 
phase carry K/K = AK of the global U(l) ^ fermion 
number charge When viewed in terms of ^, the 

transition between Hi and Hf permits a dual interpre- 
tation as a LL to band insulator quench. Correlation 
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FIG. 2: The right-moving 'supersohton.' The number den- 
sity evolution after Luttinger hquid to Mott insulator quench 
is depicted as in Fig. (TJ but here for a Gaussian initial profile 
y^Apo(a:) = Q exp (— x^/A^), with a — 0.7, now obtained 
via numerical integration of the exact bosonization result [2l[ . 
Time series for four different Q are plotted; the densities are 
normalized relative to these. The evolution is reflection sym- 
metric about X = 0. 



functions in the homogeneous quench [po{x) = 0] have 
been previously studied in Refs. (sl, Il3|. 

To characterize the post-quench dynamics, we compute 
the expectation values of the particle number (p), kinetic 
(/C) and potential {U) energy densities (the latter two 
observables are defined with respect to Hf): 



p{t,x) 

/C(t, x) 
U{t,x) 



1 

2 ^0 
i 



:vE't(t,x)vl>(t,x):|0)^^, 



2 



(8a) 

(8b) 
(8c) 



where fdg = fdg — {df)g. In these equations, 
denotes the Heisenberg picture fermion operator whose 
dynamics are generated by Hf in Eq. (j6j). U gives the 
expectation of the cosine operator in the sine Gordon 
model [Eq. (j5j)], and can be interpreted as a (squared) 
order parameter for the Mott phase. We obtain p, /C, 
and U by solving the Heisenberg equations of motion for 
^(t,x) and exploiting the bosonization map. Given an 
arbitrary initial po{x)^ we have derived exact results for p, 
/C, and U at any time t > 0, which will appear elsewhere 

m 

The exact post-quench observables in Eq. ([8|) depend 
upon po{x)^ M, and the dynamic exponent a defined via 
Eq. ([2]). The non-interacting quench with K = K has 
cr = 0, while the interacting quench {K ^ K) has <j > 0. 
We confine ourselves to the range < a < 1, for which 
the /)o(^)-dependent contributions to p, /C, and U are 
given by ultraviolet (UV) convergent integrals [2l|. At 
(7 = 1, these observables acquire logarithmic UV diver- 
gences, suggesting the onset of a sensitive dependence on 
lattice scale details. 

We now describe our main results, which concern the 
Po(^)-dependent contributions to p, /C, and U] the be- 



havior of /C and U for the homogeneous quench po = 
p{t^x) = will be discussed elsewhere [21|. The ex- 
act leading asymptotic expression for p(t, x) in the limit 
t ^ oo was already given by Eq. ([1]), above. Let us 
specialize this result to a localized initial density pro- 
file. The interacting (<j > 0) versus non-interacting 
(cr = 0) quenches yield qualitatively different behaviors. 
For the interacting quench, Eq. ([T]) implies that p{t^x) 
develops a non-dispersive response to the initial condi- 
tion for any < a < 1. For example, a Gaussian density 
bump, ^/7^Apo{x) = Qexp (— x^/A^), induces the follow- 
ing asymptotic space- time evolution for t ^ 1/M: 



p{t,x) 



2v/7rA 



e a2 



Q r(i-a) 
-2Ar(i±^) 

X e n 



{KMaff 



72 



a/2A 
V2 



7/2 



x-t' 



(9) 



where Djj{x) denotes the parabolic cylinder function, 
t' = t/K^ and we have written out the explicit form of the 
prefactor A^^ which is non- universal for a > and de- 
pends upon Ma. The naive continuum calculation gives 
Ma = 15/16. The divergence of the prefactor at <j = 1 
indicates the onset of sensitivity to the UV sector of the 
theory. 

Eq. ([9]) implies that an antecedent Gaussian density 
bump splits into right- and left-moving non-dispersive 
waves, for generic Q, A, and K ^ K {cj > ff). In the 
long time limit, the leading response is strictly linear in 
Q, with an amplitude that grows as t'^^'^ . Two Gaussian 
bumps initially separated by a distance d ^ A can be 
used to create left- and right-moving waves which pass 
through each other without changing their form [2l| . We 
dub these rigid, non-interacting density waves 'supersoli- 
tons' to distinguish them from the elementary quantum 
solitons annihilated by the fermion field ^. We have con- 
firmed the asymptotic result in Eq. (j9j) by comparing to 
numerical integration of the exact bosonization expres- 
sion for p. The supersohton is exhibited in Fig. O 

Although the precise shape of the supersohton im- 
plied by Eq. (|9|) deforms continuously with cr, it exhibits 
the same positive-negative "dipolar" peak profile for any 
< cr < 1 (see Fig. [2]). The negative density dip rep- 
resents a local evacuation of the filled Fermi sea, which 
is infinitely deep in the Luttinger model [12|. For any 
cr > 0, the integral of the second term in Eq. (|9j) over 
real x vanishes, consistent with particle number conser- 
vation. In the limit of the non-interacting quench cr ^ 0, 
the right-hand side of Eq. (|9]) vanishes; in this case, the 
response obtains entirely from subleading terms that do 
not grow with t (and conserve the particle number), but 
which we have not written here. The same is true in 
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Eq. ([T]) because 1 when a ^ 0. 

For comparison, Fig. [3] depicts the number density 
pit^x) for the case a = 0, obtained by numerical inte- 
gration of the exact result. The main message of this 
figure is that the non-interacting post-quench dynamics 
are "passive" and dispersive, depending sensitively upon 
the details of the initial inhomogeneity and showing no 
amplification phenomena. 

In the interacting quench, the supersoliton is also ob- 
served in the relative kinetic energy density, defined as 
8K[t^ X] po] = /C[t, po] — JC[t^ x; 0], shown in Fig. IH By 
contrast, we find that the potential energy density Z//(t, x) 
does not exhibit the supersoliton on top of the homoge- 
neous background it acquires after the quench. The am- 
plification in Eq. ^ does not therefore appear related to 
a Kibble-Zurek process [H^ in the order parameter. 

The physical mechanism underlying the power law in- 
homogeneity growth in Eqs. ([T]) and (j9j) can be partially 
elucidated via an analogy to the equilibrium tunneling 
density of states (TDOS) iy{(jo) in a LL [l^. Upon tun- 
neling into a one channel quantum wire at T = charac- 
terized by the Luttinger parameter the conductance 
at a bias uj = eV diminishes as iy{uj) ^ \uj\^ where a is de- 
fined as in Eq. ([2]), but with K = 1. The physics behind 
this result is a follows: The independent LL "quasiparti- 
cles" carry a fraction K of the electron charge e [17]. The 
TDOS u{uj) vanishes as ^ because a "whole" elec- 
tron must fractionalize into a large number of pieces upon 
penetrating into the LL, and this process is prohibited by 
phase space restrictions in the low bias limit. Mathemat- 
ically, the TDOS result obtains from the Fourier trans- 
form of the electron Green's function in the LL. The t^/^ 
amplification in Eq. ([T]) is rendered by a similar mecha- 
nism in the quench: an initial LL correlation function is 
convolved with an oscillatory kernel [a product of Green's 
functions resulting from the solution to the Heisenberg 



p{t,x) p{t,x) 




FIG. 3: The number density evolution as in Fig. O but for 
the non-interacting quench K = K (a = 0). The height of the 
initial bump is Q = 0.1 in the main figure and Q = 1 in the 
inset; the evolution is reflection symmetric about x = 0. Now 
there is no fractionalization of the initial LL quasiparticles 
with respect to the insulator and, consequently, the dynamics 
are simply dispersive with no supersolitons or inhomogeneity 
growth. 



51C[t,x;po] 

0.002 




FIG. 4: The kinetic energy density supersolitons. This flgure 
depicts the post-quench space-time evolution of the relative 
kinetic energy density expectation value SJC[t, x; po] (see text) 
for K = 0.77 K and an initial Gaussian density bump. 



equations of motion for The final state Hamil- 

tonian Hf introduces a scale M, by which the analog of 
the frequency u in the TDOS is the evolution interval 
M'^t. We might therefore naively expect p{t,x) ^ t^, 
with a defined by Eq. (j2j). That the leading power is cr/2 
in Eqs. ([1]) and ([9]) obtains from a cancelation of terms. 
This suggests that the immiscibility of quantum phases 
composed of quasiparticles carrying relatively fractional 
charges may underlie both the equilibrium TDOS and 
the quench amplification. 

In conclusion, we have shown that a quantum quench 
can beget a strongly inhomogeneous state, due to the 
interplay between quasiparticle fractionalization and the 
presence of a mass scale in the final state Hamiltonian. 
Fractionalization is a robust feature of ID gapless phases, 
so we expect the inhomogeneity proliferation to occur in 
many ID quantum quenches. It would be interesting 
to consider quenches to final states away from the free 
fermion LE point where (super?) soliton-soliton interac- 
tions can play a role in the dynamics. 
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